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A short review is given of the present status of the studies of genuine confinement effects in 
multiple hadron production in hard processes. 



1 Beware of Soft Confinement 



What I tell you three times is true. 

" The Hunting of the Snark" . Lewis Carroll 



Perturbative QCD (pQCD) covers orders of magnitude in the basic hard cross sections. 
Scattered clouds on the pQCD horizon (presently, a "high-pj- anomaly" and some problems 
with hadroproduction of heavy quarks and direct photons) don't define the weather. However, 
viewing the QCD landscape it is essential to remember that its "peace and quite" is deceptive. 
The closest analogy which comes to mind comparing QCD with its electroweak SM-counterpart 
is that of hell and heaven: the former being scary but entertaining, the latter — perfect but 
boring (for a detailed review see 3). 

1.1 The name of the game: MLLA 

is 1 Although this may seem a paradox, all exact sci- 
ence is dominated by the idea of approximation. 

Bertrand Russell 

In spite of the smallness of the coupling at small distances, oc s <C 1, quarks and gluons 
willingly multiply in hard interactions. This happens because the actual parameter of the 
pt expansion gets enhanced by the log of the scale of the large momentum Q applied to the 
system: a s (Q) -> a s (Q) log Q ~ 1. Such log-enhanced contributions have to be taken care of 
in all orders, giving rise to "leading- log resummations" . The structure of these contributions 
allows for a probabilistic parton interpretation. Parton multiplication in jets is described by 
the so-called Modified Leading Logarithmic Approximation (MLLA) which embodies the exact 
angular ordering resulting from coherence in multiple soft gluon radiation. As an approximation, 
MLLA is necessary for deriving asymptotically correct pt predictions. 



QCD coherence is crucial for treating particle multiplication inside jets, as well as for hadron 
flows in-between jets. Moreover, it allows the prediction of the shape of the inclusive energy 
distributions practically from the "first principles" , apart from an overall unknown normalization 
constant. IS 



1.2 Inclusive energy spectra 

It is well known that the DIS structure functions cannot be calculated perturbatively. What 
pQCD controls is the scaling violation pattern, governed by the QCD parton evolution equation 
which describes how the parton densities change with changing the scale of the transverse- 
momentum probe: 

D{x,k A _) = ^^-P{z)®D{x/z,k±_). (1) 



d In 7r 
In the Mellin moment space, f u = Jq 1 f(x)x UJ ~ 1 dx, the equation becomes algebraic, yielding 

{rk± Ah \ (x 

/ -r-TuMk))}, -y u {ct a ) = — P u ■ (2) 
Jqo k J vr 

It is the w-dependence of the input function D uj {Qq) ("initial parton distributions") that limits 
predictability of the Bjorken-x dependence of DIS cross sections. 

In the time-like jet evolution, due to the Angular Ordering, the evolution equation becomes 
non-local in the k± space: 

D(x,k±) = — P(z)®D(x/z, z-k±); [ AO: 9' = k'JzE > 8 = k±/E] . (3) 



d In k± 7r 
Using the Taylor expansion trick, 

f d 1 8 

D(x/z, z ■ k±) = exp < Inz— — — } D(x/z, k±) = z einfc -L • D , (4) 
L o In J 

the solution in the moment space comes out similar to that for the DIS case, Eq. [| but now 
with an operator as an argument of the splitting function P: 

-=55!-- 



This leads to the differential equation 
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Since we are interested in the small-x region, the essential moments are small, uCl. For the 
sake of illustration let us keep only the most singular piece of the splitting function (DLA) , 

Then Eq. |] immediately gives a quadratic equation for the anomalous dimension, 



2N c a s m ( a 2 „\ 
( w + 7 ) 7 C -JL + Q M = 0. (8) 



NB: It suffices to use the next-to- leading approximation to the splitting function, P^ ~ 2N c /uj—a, 
with a = ll-/V c /6 + nj/(3A r ( ?), and to keep the leading correction coming from differentiation of 
the running coupling in Eqs. [f| [8|, to get the more accurate MLLA anomalous dimension 7 W . 



The leading anomalous dimension following from Eq. B is 



7 = | (-1 + v/l + SiW^ 2 ) • (9) 

When expanded to the first order in a s , it coincides with that for the space- like evolution, 
7 ~ OLgj'K ■ Pu, with P given in Eq. [7[ 

The time-like DLA anomalous dimension Eq. |9| (as well as its MLLA improved version) has 
a curious property. Namely, in a sharp contrast with the DIS case, it allows the momentum 
integral in Eq. ^ to be extended to very small scales. Even integrating down to Qq = A, the 
position of the "Landau pole" in the coupling, one gets a finite answer for the distribution (the 
so-called limiting spectrum), simply because the y/a s (k) singularity happens to be integrable! 

It would have been a bad taste to actually trust this formal integrability, since the very 
perturbative approach to the problem (selection of dominant contributions, parton evolution 
picture, etc) relied on a s being a numerically small parameter. However, the important thing 
is that, due to time-like coherence effects, the (still perturbative but "smallish") scales, where 
a s (k) u 2 , contribute to 7 basically in a co-independent way, 7 + to/2 oc y/a s (k) ^ f(w)- This 
means that "smallish" momentum scales k affect only an overall normalization without affecting 
the shape of the ^-distribution. Since such is the role of the "smallish" scales, it is natural to 
expect the same for the truly small — non-perturbative — scales where the partons transform 
into the final hadxons. This idea has been formulated as a hypothesis of local parton-hadron 
duality (LPHD)H 

According to LPHD, the x-shape of the so-called "limiting" spectrum which one obtains by 
formally setting Qq = A in the parton evolution equations, should be mathematically similar 
to that of the inclusive hadron distribution. Another essential property is that the "conversion 
coefficient" should be a true constant independent of the hardness of the process producing the 
jet under consideration. Starting from the LEP-I epoch, this "prediction" stood up to scrutiny 
by e + e~, DIS and Tevatron experiments. 
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Fig.l: CDF hump-backed plateau versus an ana- 
lytic MLLA prediction for the yield of secondary 
partons (soft gluons).0 
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Fig. 2: Position of the maximum in the in- 
clusive energy spectra versus a parameter- 
free MLLA prediction. § 



The message is, that "brave gluon counting" , that is applying the pQCD language all the way 
down to very small transverse momentum scales, indeed reproduces the x- and Q-dependence 
of the observed inclusive energy spectra of charged hadrons (pions) in jets. 

Even such a tiny (subleading) effect as an envisaged difference in the position of the maxima 
in quark- and gluon-initiated humps seems to have been verified, 15 years later, by the recent 
DELPHI analysis §§ 

1.3 Inter- jet particle flows 

"Can you do addition?" the White Queen asked. "What's one 
and one and one and one and one and one and one and one 
and one and one?" "I don't know," said Alice. "I lost count." 

Through the Looking Glass 

Even more striking is miraculously successful role of gluons in predicting the pattern of 
hadron multiplicity flows in the inter-jet regions — realm of various string j drag effects. It isn't 
strange at all that with gluons one can get, e.g., 1 + 1 = 2 while 1 + 1 + | = j^, which is a 
simple radiophysics of composite antennas, or quantum mechanics of conserved colour charges. 



This particular example of "quantum arith- 
metics" has to do with comparison of hadron 
flows in the inter-quark valleys in qqj and qqg 
(3-jet) events. The first equation describes the 
density of soft gluon radiation produced by two 
quarks in a qq'j event, with 1 standing for the 
colour quark charge. 

Replacing the colour-blind photon by a gluon 
one gets an additional emitter with the relative 
strength 9/4, as shown in the l.h.s. of the second 
equation. The resulting soft gluon yield in the 
qq direction, however, decreases substantially as 
a result of destructive interference between three 
elements of a composite colour antenna. 
Nothing particularly strange, you might say. 
What is rather strange, though, is that this 
naive perturbative wisdom is being impressed 
upon junky 100-200 MeV pions which dominate 
hadron flows between jets in the present-day ex- 
periments such as the OPAL study shown in 
Fig. 3. 

Another amazing test of this sort was provided by the DELPHI measurement of the multiplicity 
of "(low energetic) tracks emitted perpendicular to the eveny plane" in 3-jet events^ which has 
been found to obediently follow a simple pt prediction based on coherent soft gluon radiationQ 

• The colour field that an ensemble of hard primary partons (parton antenna) develops, 
determines, on the one-to-one basis, the structure of final flows of hadrons. 

• The Poynting vector of the colour field gets translated into the hadron pointing vector 
without any visible reshuffling of particle momenta at the "hadronisation stage" . 

When viewed globally, confinement is about renaming a flying-away quark into a flying-away 
pion rather than about forces pulling quarks together. 
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Fig. 3: Comparison of particle flows in the 
qq valley in qq'j and qqg 3-jet events versus 
a parameter-free analytic prediction based 
on the soft gluon radiation pattern, i 



1.4 Gluons and Gluers 



Definition: a Gluer is a miserable gluon which hasn't got enough time to truly behave like one 
because its hadronization time is comparable with its formation time, iform. — w/k"j_ ~ ihadr. — 
ujR 2 onf . Contrary to respectful pt gluons born with small transverse size, k± S> -R^nf , gluers 
are not "partons": they do not participate in perturbative cascading (don't multiply). Accord- 
ing to the above definition, gluers have finite transverse momenta (though may have arbitrarily 
large energies). Having transverse momenta of the order of inverse confinement scale puts gluers 
on the borderline of applicability of pt language, since their interaction strength is potentially 
large, cn s (R~^ nf ) ~ 1. Role of gluers is to provide comfortable conditions for blanching colour 
parton ensembles (jets) produced in hard interactions, locally in the configuration space. Gluer 
formation is a signal of hadronization process taking place in a given space-time region. A label 
to put on the gluer concept might be — 11 A gluer formed ~ a hadron born" . An Idea emerges: 
To relate (uncalculable) Non-Perturbative corrections to (calculable) Perturbative cross sec- 
tions/observables with intensity of gluer emission (a s in the infrared domain). 



2 Power Games 99 



PT-calculable observables are Collinear-and-InfraRed-Safe (CIS) observables, those which can 
be calculated in terms of quarks and gluons without encountering either collinear (zero-mass 
quark, gluon) or soft (gluon) divergences. Gluers' contributions to such observables are sup- 
pressed and are being rightfully neglected in the pure pt ( "logarithmic" ) approximation. These 
contributions are inversly proportional to a certain power of the hardness scale (modulo logs), 
(5cr NP ja oc log 9 Q/Q 2p . The corresponding observable-dependent exponents can be inferred from 
the analysis of an intrinsic uncertainty in summing up the pt series (infra-red renormalons, for 
an extensive review seet3). 

Adopting the concept of universality of np phenomena one can predict the ratio of the 
magnitudes of power corrections to different observables belonging to the same {p, q} class. 

The PT-approach exploiting gluers allows to go one step further, namely to relate absolute 
magnitudes of genuine np contributions to CIS observables with the intensity of gluer radiation, 
i.e. the "QCD coupling" at small transverse momentum scales. 



2.1 Phenomenology 



For example, DIS structure functions are 
expected to deviate from their perturbative 
Q 2 dependence by terms generally behaving 
like l/Q 2 ("twist 4"): 



F 2 (x, Q^) ~ F2 (x, Q ) [l + D 2 (x,Q 2 )/Q 2 

Comparison of the Power Game predic- 
tion0 with the data@ allows one to extract 
the value of the characteristic np parameter 
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Fig.4: x-dependence of the l/Q 2 contribu- 
tion to F2 

A2 being fixed, a parameter-free prediction emerges then for the l/Q 2 suppressed contribution 
to F3 shown by the dashed curve. 

Another example is provided by a variety of observables including jet shapes (Thrust, C- 
parameter, jet Broadenings, Oblateness), energy-energy correlation (EEC), ctl, etc. which belong 



to the p = 1/2 class and thus exhibit numerically large np deviations. 
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Fig. 5: Mean Thrust and C-parameter in e + e~ annihilation El 
A pure phenomenological study of the deviation of the mean Thrust and C-parameter values 
from the corresponding two-loop pt predictions shown by dotted lines in Fig. 5 hints at 

ci(l-r) (NP) = (l-T) - (1-T) (PT) ~ lGeV/Q, 

5(C) (NP) = (C) - (C) (PT) ~ 4GeV/Q, 

with the power game bet being instead 5 (C) ( - NP ^/(5 (1 — T)^ NpS> = 3tt/2. 

2.2 Universality of confinement effects in jet shapes 



The Power Game grew muscles when 
it was realised that it can be played 
not only with the Q-dependence of the 
means at stake. The distributions of 
shape variables were shown to be 
subject to a 1/Q shift, by that very 
amount that describes the genuine np 
contribution to the mean value of the 
corresponding jet shape variable. 
For example, the C-parameter distri- 
bution (for the values of C not too 
close to zero) can be obtained by 
simply shifting the corresponding all- 
order-resummed purely perturbative 
spectrum by an amount inverse pro- 
portional to Q, 



1/a do/dC 



(preliminary) 



1 da 

adC 



(C)c 



Q J 



The corresponding result of a recent 
JADE analysis is shown in Fig. 6. 
The same shift prescription, and sim- 
ilar high quality description, hold for 
other CIS jet observables like Thrust. 
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Fig. 6: C-parameter distributionsL£l versus PT-spectra 
shifted by 5C ~ 4 GeV/Q 



A thrilling story of one important exception was told in Vancouverll^l Jet Broadening (s) 
denned as a sum of the moduli of transverse momenta of particles in jet(s) (wrt the Thrust 
axis) was first predicted to have a log Q-enhanced np shift, since this np contribution to B was 
naturally thought to accumulate gluers with ramdities up to log Q. 

The data however simply could not stand itl 17 H 18 l Fits based on the log Q-enhanced shift were 
bad and produced too small a value of a s (Mz), and the np parameter «o inconsistent with that 
extracted from analyses of the Thrust and C-parameter means and distributions. 

Tragic consequences for the universality belief seemed imminent. 



2.3 Broadening: tragedy, catharsis, lessons 

vs 1 Raffiniert ist der Herrgott, aber nicht boesartig 

A. Einstein 

Catharsis came with recognition of the fact that the Broadening measure (B) is more sen- 
sitive to quasi-collinear emissions than other jet shapes, and is therefore strongly affected by 
an interplay between pt and np radiation effects. With account of the omnipresent pt gluon 
radiation, the direction of the quark that forms the jet under consideration can no longer be 
equated with the direction of the Thrust axis (employed in the definition of B). As a result of 
this interplay, the hadron distribution was found to be not only shifted but also squeezed with 
respect to its pt counterpart. 




Fig. 7: Perturbative (dashed) and NP-shifted/squeezed Total Broadening distributions 
Three lessons can be drawn from the Broadening drama. 



Pedagogical lesson the Broadenings taught, was that of the importance of keeping an eye on pt gluons 
when discussing effects of np gluers. An example of a powerful interplay between the two 
sectors was recently given by the study of the energy-energy correlation in e + e _ in the 
back-to-back kinematics^ The leading 1/Q np contribution was shown to be promoted by 
pt radiation effects to a much slower falling correction, Q-°-32-o.36_ 

Physical output of the proper theoretical treatment was restoration of the universality picture: 
within a reasonable 20% margin, the np parameters extracted from T, C and B means 
and distributions were found to be the same. 

Gnostic output was also encouraging. Phenomenology of np contributions to jet shapes has shown 
that it is a robust field with a high discriminative power: it does not allow one to be misled 
by theorists. 



2.4 A step forward: "shape functions" 

A strong push is being given to the Power Game by the notion of Shape Function(s) introduced 
by G. Korchemsky and G. StermanEil 




Introducing the distribution describing the 
power shift on an event by event basis (shape 
function) makes it possible to lift off the condi- 
tion 1 — T ^> (1 — T)^ P \ A simple physically 
motivated ansatz for such a distribution for the 
Thrust case, 



1 da(l-T) 



O~tot 

with 
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Q(i-T) da 
def(e) — 



1-T 



dT 



2(e/A) a ~ 1 
AT(f) 



exp 




produces a remarkable fit to hadron data shown 
in Fig. 8. 



Fig.8: Thrust distributions at Q=U, 22, 
35, 44, 55, 91, 133 and 161 GeV versus 
Power Game prediction Ell 

Shape functions for different jet shapes can be related with certain characteristics of the energy- 
momentum flow at the hadronization stage, specific for a given observable. 



2.5 Universality problem 

A detailed discussion of the main problems one faces inestablishing the rules of the Power 
Game can be found in the proceedings of the 1998 ICHEPt^l These problems include separation 
of power corrections coming from the infrared region from those determined by the ultraviolet 
physics, merging (in a renormalon-free manner) the PT and NP contributions to the full answer, 
the problem of splitting the magnitude of the power term into an observable-dependent pt- 
calculable factor and a universal np parameter. 

The key question is whether the latter is really universal. The whole game would have little 
sense if it were not. Allowing each observable to have a private fitting parameter we would not 
learn much about the way confinement acts in hadronizing ensembles of partons produced in 
hard interactions. 

Reasonable doubt was expressed in a seminal paper by Nason and Seymouril as to whether 
universality can be expected to hold for jet shapes which are not truly inclusive observables. 
The configuration of offspring partons in the gluon decay matters for jet shapes, so that the 
value of the power term may be affected, in an observable-dependent way, beyond the leading 
level in a s (which a priori is no longer a small parameter since the characteristic momentum 
scale is low). 

Analyses of two-loop effects in 1/Q suppressed contributions have been carried out for jet 
shapes in e + e~ annihilation and DIS. The output proved to be surprisingly simple. It was shown 
that there exists a definite prescription for defining the so-called "naive" one-loop estimate of the 
magnitude of the power contribution, such that the two-loop effects of non-inclusiveness of jet 
shapes reduce to a universal, observable-independent, renormalisation of the "naive" answer by 
the number known as the "Milan factor" B3'Ej This is true for the np contributions in the thrust, 
invariant jet mass, C-parameter and broadening distributions, for the energy-energy correlation 
measure, as well as for other observables subject to linear in 1/Q confinement effects. 

It is probably the striking simplicity of the resulting prescription to be blamed for apparently 
cold reception the "Milan factor" enjoyed among theoreticians. 

Verification of the Milan factor prescription is underway. M. Dasgupta, L. Magnea and 
G. Smye have undertaken the project of explicitly calculating the two-loop effects in the np 



contribution to the C-parameter distributionCJ The analytical result they are coming up with 
has verified the key simplification used in the original derivation of the Milan factor namely, the 
soft gluon approximation. This is good news. The not-so-good news is that the final expressions 
for M do differ... 

3 Milan factor 2000 

The Power Game as a new theoretical instrument emerged from its toddler years but has not 
yet reached respectable teens. It is understandable that, being both predictive and verifiable 
(the qualities almost extinct nowadays), it attracted a lot of attention and was developing, in 
its early days, on a week-to-week (if not a day-to-day) basis. Accelerated childhood tends to be 
marked by bruises, on the child's part, and by troubles on the parents'. 

A partial history of misconceptions the advocates of the Power Game had to muddle through 
can be found inE3. Now we are in a position to enrich this history with a miscalculation. An 
unfortunate omission of a trivial factor in the two-parton phase space resulted in a wrong value 
originally derived for A4: the so-called "non- inclusive" contribution to the Milan factor, r m , has 
to be multiplied by a factor of 2. As a result, 

M = l + r i „ + r ni ^l + 3.299C A //3o + 2 x(-0.862C A -0.052n / )//3o 

= 1 + (1.575CU - O.KMn/V/Jo = 1.49 for n f = 3 , (instead of = 1.8) . (10) 

The nj-part of the corrected Milan factor Eq. [l(] agrees withil and, as the authors point out, 
also solves the longstanding discrepancy with the explicit two-loop calculation of the "Abelian" 
(^/-dependent) correction to singlet e + e~fragmentation functions (cr^) which was carried out 
by M. Beneke, V. Braun, and L. Magneac3 

Refitting jet shape data with the corrected A4 lies ahead. It will drive up the np parameter 
ao by about 10% but will change neither a s nor the present status of the universality pattern. 

The situation with universality these days can be viewed as satisfactory. It is far from 
perfect, however. In particular, there seems to be a conceptual problem with describing the 
means and distributions of those specific jet variables that deal with a certain single jet rather 
than the event as a whole. The known cases this remark applies to, are the Heavy jet mass 
and the Wide jet broadening. An adequate game strategy for dealing with such (less inclusive) 
observables remains to be found. 

A last remark is due concerning the title "Power Games" . An ideology and technologies are 
being developed for describing genuine confinement effects in various global characteristics of 
multi-particle production. I believe there was a good reason for calling it a "game". To really 
enjoy playing one has to follow the rules (which, by the way, does not contradict the fact that 
some entertaining games intrinsically embody bluff). 

In the present context, "the rules" means equating 
"pt" with the two-loop prediction and looking upon 
the rest as being "np". The boundary between pt 
and np physics is, to a large extent, a matter of con- 
vention. In particular, includingan additional loop 
into a "pt prediction" (see, e.g. E^) or redefining it, 
say, with use of the Borel wisdompj inevitably affects 
the magnitude of a "genuine np contribution". Such 0.1 
an elusive behaviour of np effects may appear espe- 
cially confusing in jet shape phenomenology where, 

according to the Sterman's lemma^l the NNLO (a|) q [Gevj 

effects are perfectly capable of mimicking the l/Q Illustration of the Sterman's lemma 
behaviour. 
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